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We derive an uncertainty relation for the energy density and pressure of a 
quantum scalar field in a time- dependent, homogeneous and isotropic, classical 
background, which implies the existence of large fluctuations comparable to 
their vacuum expectation values. A similar uncertainty relation is known to 
hold for the field square since the field can be viewed as a Gaussian random 
variable. We discuss possible implications of these results for the reheating 
process in scalar field driven inflationary models, where reheating is achieved 
by the decay of the coherently oscillating infiaton field. Specifically we argue 
that the evolution after backreaction can seriously be altered by the existence 
of these fluctuations. For example, in one model the coherence of the infiaton 
oscillations is found to be completely lost in a very short time after backreaction 
starts. Therefore we argue that entering a smooth phase in thermal equilibrium 
is questionable in such models and reheating might destroy the smoothness 
attained by inflation. 
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Our universe is very special. The smoothness of the cosmic microwave background (CMB) radiation 
is the prominent evidence for this assertion. It is remarkable that CMB photons were in thermal 
equilibrium on scales much larger than the horizon size at the time of decoupling. Therefore, causality 
precludes thermalization of the CMB radiation by local interactions and the problem we are facing 
is completely different than a box of gas reaching thermal equilibrium. Actually, there is a deeper 
difficulty since in the standard cosmological model without inflation, most of the CMB photons did 
not have a chance to causally communicate at all. 

Inflation claims to solve this problem by exponentially enlarging a small causal patch, which 
later becomes our observed universe. However, when inflation ends one does not observe thermal 
equilibrium right away. Rather, in the scalar field driven models for example, one finds an extremely 
smooth, flat universe flUed only with a coherently oscillating inflaton fleld; all rest that existed before 
inflation are red-shifted away by the exponential expansion. Thermal equilibrium must be reached 
by the interactions accompanying the decay of the inflaton fleld. 

This shows that one still needs to understand the thermalization process in an expanding universe, 
which is a very non-trivial and difficult problem to study (see [1]). To have a correct picture in this 
setup, it is very crucial to determine the state of the universe just before the thermalization starts. 
Our aim in this essay is to point out some important features of quantum particle creation in a 
classical background, which must be taken into account during reheating and before thermalization. 

The particle creation process during reheating is usually studied using Fourier decomposition and 
momentum modes (see e.g. [2]). While it is perfectly legitimate to use momentum modes, or in 
essence any complete set of modes, one must be careful about a few potential issues. Firstly, causahty 
requires the particle creation process to take place independently in each horizon. Thus, in physically 
interpreting the creation of a globally defined momentum mode one must ensure that locality is 
not broken down by hand. Secondly, there is a subtlety in the definition of the number density 
of a completely dislocalized momentum mode. Finally, to make sense of local physical quantities 
a suitable regularization must be utilized and this may also change some results obtained in the 
momentum space. An important example of this kind is the modification of the power spectrum 
by adiabatic regularization, studied recently in [3]. All these issues can be bypassed if the particle 
creation effects are described by giving the vacuum expectation value of the energy-momentum 
tensor, which we adapt in the following. 

Consider now a real scalar field x propagating in a Friedmann-Robertson- Walker metric 

ds^ ^ -dt^ + a^{dx^ + dy^ + dz^), (1) 

which has the following action 

'5 = -^/v^[(Vxf + MV]- (2) 

We assume that the mass parameter M may also depend on time: M = M{t), and therefore particle 
creation effects occur due to the time dependence of the scale factor a and the externally varying 
mass parameter M. For quantization, it is convenient to define X — a^/^x? introduce the Fourier 
modes and time-independent creation- annihilation operators as 



akXfce-^'^-^ + 4x*e^^-^' , (3) 



where [ok, a^,] — (5(k — k'), Xk obeys the Wronskian condition XkX^ — X^Xk — i and Xk+^^l-^k — 
with the frequency defined as ul = + ^ ~ f-^^^ — f-^- 

The ground state of the system is defined as ak|0 >= 0. It is a straightforward exercise to calculate 
the vacuum expectation value of the energy momentum tensor < T^i, >=< 0|T^,^|0 > corresponding 
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to the x-particles created on the background. By noting that the energy density and the pressure 
of the x-particles are given by 

P{t: x) = - l9''{^^x){^JX) - \m\\ (4) 



one can obtain 



<p>=T + l^ + G, <P >=T -V -\g, (5) 



where T, V and G are given by 

2(27ra)3y ' 2 ' 



M2 

y = 7T77^ / \Xk\\ (6) 
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Also, it is of some interest to determine the field variance < x^ >—< 0|x^(t,x)|0 >, which can be 
found as 



It is clear that T and V can be interpreted as the kinetic and the potential energies, respectively, 
and G can be viewed as the energy stored in the gradient of the x-ficld. In general all the above 
expectation values diverge and thus a suitable regularization must be applied. After regularization, 
however, the functions T, V , G and < > can be made finite. Note that since the background is 
homogenous the vacuum expectation values (5) and (7) do not depend on the spatial coordinates. 

The above expressions only give the average values of the corresponding physical quantities and 
one would expect to determine (quantum) fluctuations about these averages. In [4], we explicitly 
calculate them. For example, after defining the fluctuation operators 5p — p(i,x)— < p > and 
5P — P{t,x)— < P >, a relatively long but straightforward calculation, which uses (3) and (4), 
gives 

20 

<5p'^> + < 8P'^ >= 4T^ + 4^2 + — (8) 

27 

From (5) this result shows that the deviations in the energy density and the pressure are always 
of the order of the average energy density < p >. Similarly, although their vacuum expectation 
values vanish one can see that the momentum density and the stress tensor have non-zero variances 
proportional to TG and G^, respectively [4]. Thus, if T and G are not small compared to V, 
the deviations in the momentum density and the stress can also be as large as < p >. Eq. (8) 
can be viewed as an uncertainty relation involving energy density and pressure, where the order of 
magnitude of the uncertainty is fixed by the average energy density. 

It is easy to see that there also exists order one fluctuations in x^, i-e. the field variance has a 
large variance. Defining 5x^ = x^(^)X;)— < x^ >, one can easily calculate 



< {Sx'r > = V2<x'>. (9) 
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This result directly follows from the fact that x) is a Gaussian random variable defined at the 
point X and (9) is true for any Gaussian distribution. 

The fluctuations in (8) and (9), which are actually defined pointwise, should be interpreted as the 
statistical averages over different points in space at a given time. However, it is crucial to realize 
that the field variables at nearby points are not statistically independent. The (comoving) size of 
a typical region containing correlated field variables is given by the correlation length ,^c, which 
is the minimum value of |x| such that the two-point function vanishes < x(^)X)x(t,0) >= 0. The 
correlation length is one of the most important parameters characterizing quantum fluctuations. One 
can imagine the space to be divided into regions of typical size ^c-, and in each region the physical 
quantities p, P and can be thought to be uniformly distributed. The variations given in (8) and 
(9) characterize changes from one such region to another. 

It is interesting to compare these findings with [5, 6] where we study the particle creation process 
using a complete orthonormal family of localized wave-packets and find out existence of fiuctuations 
in the number and the energy densities of particles produced in any given volume V . As noted in 
[5, 6], if the volume V is sufficiently large, the relative deviations of these perturbations fall down 
like 1 / \/V . This behavior can easily be understood since the total deviation is due to V/ vq number 
of statistically independent random variables, where Vq is the volume occupied by each random 
variable, which can be fixed in terms of the correlation length as Vq — 

One may now wonder if the presence of these fiuctuations can have any cosmological significance. 
Here, we consider the reheating process in a single scalar field model. Assume that near its minimum 
the potential of the inflaton field takes the form ^rn^cj?^ where m is the infiaton mass. In that case 
the inflaton and the scale factor of the universe a{t) can be determined as 

= $(t)sin(mi), ~ — , a^[ — ] . (10) 

t \to/ 

The pressureless dust equivalent expansion in (10) can be understood as follows: As the inflaton 
oscillates about its minimum, the kinetic and the potential energy terms consecutively dominate the 
energy-momentum tensor. While the kinetic energy dominance can be described by the equation of 

state P = p the potential energy dominance implies P = —p. Thus on the average one gets P = 0. 

We further focus on a model where the inflaton decays to a minimally coupled x-boson via the 
interaction 

1 1 

Ant = -^g'^(f>'^X^ or Cint = -^cr(l>X^- (11) 

The first term can give a decay in the broad parametric resonance regime and the second one 
is usually considered in the perturbative decay. In either case, at least in the first stage until 
backreaction becomes important the decay process can be described as the x-^artic/e creation in a 
time dependent classical background. One can thus apply the above formulas about particle creation 
where the time-dependent mass parameter M can be fixed as = g'^cp'^ or = a(j). 

In [4] , we determine the correlation length in these two different reheating models corresponding 
to the decay in the broad parametric resonance regime and in the perturbation theory. In the former 
case, is fixed by the (comoving) momentum scale /c* of the first instability band as ~ ^/k*- On 
the other hand, in perturbation theory is determined up to the scalings due to the expansion of the 
universe as ~ V"^- obtaining both results one simply uses the spectrum of particles produced 
during the decay in calculating the two-point function. Importantly, in each case the correlation 
length turns out be smaller than the Hubble scale. 

The above formulas can be used until backreaction effects become important. Prior to backreaction, 
linearity is essential and the momentum modes evolve independently without disturbing locality. 
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However, as the created particles start influencing the background, the hnearity is lost and one 
can no longer treat momentum modes separately. Namely, to determine how the created particles 
alter the evolution of the background, which necessarily obeys local field equations, the efi^ects of 
all the modes must be summed up. Therefore, it is crucial to determine how physical quantities 
corresponding to particle creation behave in the real conflguration space. 

Consider, for example, the model given by the first interaction term in (11), which modifies the 
inflaton field equation as 

+ 3H^ + m^ct) - g'^{didj(f)) + c/^xV = 0, (12) 

where wc also introduce the gradient term for completeness. The backreaction effects become impor- 
tant when the ^-field grows so that g^x^ ~ the Hartree approximation one replaces term 
in (12) with the expectation value < >■ Since < > does not depend on spatial coordinates 
and since initially (i.e. just after the end of inflation) the 0-fleld is homogenous, only the zero mode 
continues to exist. Thus, in this approximation one can ignore the spatial derivatives in (12) and 
the 0-fleld remains to be homogenous. 

However, as noted above in (9) there exists large fluctuations in x^, comparable to its vacuum 
expectation value < x^ > ^li^s the " actual value" of x^ varies appreciably on scales larger than 
the correlation length ^c- To determine the evolution of the inflaton fleld in the presence of these 
fluctuations, one can examine (12) in regions of volume , in which x^ can be taken as uniform, 
and try to glue these local results suitably. It is clear that if only the zero mode survives even when 
the fluctuations are taken into account, then this should be a good approximation. 

Prom (12), the frequency of the oscillations in the i'th region is given w^j) = + g'^xh where Xi 
denotes the value of x^ in that region and we ignore the expansion of the universe since in general 
m ^ H during reheating. This shows that due to fluctuations in x^ given by (9) the frequencies 
also start to differ as 

5uj 5x^ , 
u x^ 

In a very short time nearly corresponding to a single average oscillation t ~ l/u ~ 1/m, the 
oscillations of the inflaton fleld in different regions become completely out of phase. In that case, 
it is no longer permissible to ignore spatial derivatives in (12) and the whole reheating dynamics 
changes completely. 

Actually, the story is more complicated since x^ is a random variable. Consider, for instance, the 
situation when g^ < X^ > is about 5% percent of vr? or so, and thus backreaction effects can be 
ignored on the average. However, even at that time one can find regions in which g^')^ ~ w? and 
thus coherence of the oscillations is lost near such regions. In other words, one cannot actually talk 
about a definite time after which backreaction effects become important. 

The above conclusion is specific to the interaction g^'^'^y^ and one may wonder if it also holds for 
other models in general. For example, when the infiaton decays through the interaction cr0x^) 
backreaction term does not modify the frequency of the inflaton oscillations. But this time, there 
appears a non-homogenous source term ax^ in the inflaton field equation which alters the amplitude 
of the oscillations locally. One can see that the fluctuations in x^ induce order one changes in the 
amplitude and again the gradient terms cannot be neglected in the inflaton fleld equation. 

Thus, we flnd out that in reheating models characterized by the interactions (11), the smooth- 
ness of the inflaton background is completely lost and the subsequent evolution must be determined 
by taking into account the field derivatives. One may think that since the length scale of inho- 
mogeneities is small compared to the cosmologically relevant scales today, they cannot have any 
cosmological imprints. However, it is not possible to view these fiuctuations as small perturbations 



(13) 
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on a homogenous background. Namely, they are expected to alter the evolution of the background in 
a very non-trivial way (for instance the gradient terms may start dominating the energy-momentum 
tensor) and the universe may transform into a state which is completely different than a homogenous 
and isotropic Robertson- Walker space. 

One may wonder what can be said for other reheating models in general. Until backreaction 
effects become significant, (8) holds in any model of reheating provided that the decay process can 
be modeled as the quantum particle creation in a classical background. Therefore, in such models 
there inevitably exists order one fluctuations in the energy density and the pressure. As is known, in 
some cases the universe can still be described as a perturbed Robertson- Walker metric even though 
5p/p > 0{1). However, the fluctuations corresponding to the particle creation have interesting 
properties. For example, it is not possible to characterize them by an equation of state [4] and 
consequently the combination p + 3P, which determine a, would change locally causing different 
regions to accelerate with different rates. Moreover, since the size of each region is subhorizon 
and there exists large pressure gradients (as shown in [4]), the gravitational collapse can occur 
exponentially fast as in the case of Jeans instability in Newtonian theory. As a result, one expects 
these fluctuations to grow and the universe may not enter into a smooth phase in thermal equilibrium. 

In the above discussion we left over some important issues which must actually be studied in real- 
istic scenarios. It is known that for completeness the metric and the inflaton fluctuations must also 
be incorporated in the linearized fleld equations which may alter the particle creation characteristics 
[7]. In that case, however, it is not surprising to see the existence of large fluctuations of these flelds 
as well. On the other hand, one may also need to consider the rescattering or the thermalization 
effects, which may work for uniformization. There are also some conceptual problems one should 
solve. For example, in (12) while </> is a classical fleld it is not obvious what should stand for. 
Here, we take the view that it should be replaced by a classical realization which would agree with 
quantum expectation values. However, especially after backreaction sets in, it is not very clear how 
to determine the evolution of the fields since quantum and classical objects start mixing. All these 
issues must be resolved before concluding that reheating ends with a perfectly smooth universe in 
thermal equilibrium. On the contrary, our findings indicate that the opposite is likely to occur. 
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